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Abstract 

A controversy is raging concerning the most physically meaningful definition of quark and gluon 
angular momentum, a question of importance in understanding the internal structure of the nu- 
cleon. There are, in my, opinion, two principal candidates: the Belinfante version, JbeU which is 
gauge invariant, and the Canonical version, Jcani which is non-gauge invariant, but which has a 
physically intuitive structure. I demonstrate the wholly unexpected result that for quarks in a spin 
1/2 nucleon the expectation value of J^i (quark) is equal to the expectation value of J can (quark) 
taken in the light-cone gauge A + = 0, with an analogous result for gluons. 

PACS numbers: 11.15.-q, 12.20.-m,12.38.Aw,12.38.Bx, 13.88. +c,13.60.Hb,14.20.Dh 
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Perhaps somewhat surprisingly, it has turned out that defining what is meant by quark 
and gluon angular momentum is a far from trivial matter, and there has been a major 
controversy about this, culminating in a Workshop at the INT in Seattle in February 2012 



for access to the literature see http://www.int.washington.edu/PROGRAMS/12-49w and 



1] . One of the issues, the splitting of the gluon (or photon) angular momentum into orbital 
and spin parts, long regarded as impossible in a gauge invariant way, is highlighted by the 
fact that we nonetheless measure AG(x), the gluon polarization in the nucleon. Several 
suggestions for achieving a gauge invariant split have been made {2 but all involve the 
introduction of non-local fields. I believe that there are really only two fundamental variants 
of the angular momentum operators: the Canonical version which emerges directly from 
Noether's theorem, but which is not gauge invariant, and the Belinfante version, which is 
obtained from the Canonical density by adding a spatial divergence, and which is gauge 
invariant. 

It is well known that for the total angular momentum of a system the expectation values of 
Jbei and J C an are equal in a spatially localized state, since the operators differ by an integral 
of a spatial divergence. The situation is quite different when the total angular momentum is 
split into contributions from separate constituents, e.g. quarks and gluons. Thus for quarks, 
Jbei, q and J C an,q do not differ by an integral of a spatial divergence and there is no a priori 
reason to suppose that their expectation values will be equal; and analogously for the gluon 
contribution. 

I shall show what is therefore a remarkable and unexpected result, that in a spinl/2 nucleon 
the expectation value of the quark angular momentum Jbei,q is equal to the expectation 
value of J ca n,q evaluated in the light-cone gauge A + = 0, and similarly for the gluon angular 
momentum. 

The three angular momentum operators J 1 are related to the angular momentum density 
Mw(x) by 

j = \tijkJ jk J jk = j <? x M° jk (t, x )- (!) 

and are independent of time for the total, conserved angular momentum. The Belinfante 
angular momentum density is given in terms of the energy momentum density by 

MZ°{x)=xOTZ-x°TZ. (2) 
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Substituting Eq. fl2]) into Eq. (pjj and using translational invariance : 
J d 3 x(P|^(x)-^(x)|P) 

= J d s xx* (P\e ip - x T^(0)e- ip - x \P)-(i^j) 

= J d 3 x^(P|T°i(0)|P)-(z^j)- (3) 

Since the matrix element in Eq. ([3]) is independent of x the integral is totally ambiguous, 
being either infinite or, by symmetry, zero. There are two ways to deal with this: using a 
careful wave-packet analysis as was done by Bakker, Leader and Trueman (BLT) [5[ or using 
a trick due to Jaffe and Manohar (JM) j(J. In this paper we shall follow the latter approach. 
Consider 

M^ x (P,k,P',S,S') = J d 4 xe ik - x (P',S'\M^ x (x)\P, S) (4) 

which is designed to regulate the ambiguity in Eq. fl3]). Each state is specified by its 4- 
momentum P and a covariant spin vector S, normalized to S 2 = —M 2 . S and S' are 
functions of the different momenta P and P' respectively, but of the same rest-frame spin 
vector. One may wonder why a 4-dimensional Fourier Transform is introduced in dealing 
with a 3-dimensional integral. The reason is that M.^ uX (P, k, P', S, S') appears to transform 
as a covariant Lorentz tensor, but that is an illusion, because the non-forward matrix element 
of a tensor operator is not a tensor. This misunderstanding is partly responsible for an error 
in j^l (for a detailed explanation and consequences of the error, see 5|). Note also that 
for the validity of later manipulations P' should be treated as an independent variable and 
should not be put equal to P + k. Eventually the limits fc M — > and P' — > P are taken, and 
this will yield an expression for the matrix elements of J jfc up to a trivial factor. Firstly 

lim p hm p M 0]k (P, k, P', 5, S') = j dte lkot J d 3 x{P', S'\M^ k (x)\P, S) 

= 2n5(k )(P,S\Jt k \P,S) (5) 

and then finally 

lim lim M Ojk (P,k,P',S,S') = 27r8(0)(P,S\Ji k \P,S) (6) 

and the innocuous factor 2tt5(0) cancels out later in the calculation. In summary the JM 
trick is very useful provided one is careful not to give M 0l/X (P, k,P',S,S') incorrect trans- 
formation properties. As shown by the wave-packet analysis of BLT, Ai° 3k (P,k,P',S,S') 



contains two terms, one referring to the internal angular momentum of the nucleon, the 
other to the angular momentum of the packet about the origin of coordinates. Only the first 
term is of interest and after some manipulation yields 



(P,S\J&\P,S) 
i(2vr)¥ 3 (0) 







P + A/2,S\T™\P-A/2,S)-(j^k) 



8Aa 



(7) 



Clearly, from Eq. ((Tj) we only need an expansion of ( P + A/2, S\ \ P — A/2, S ) to first 



order in A, but this is a little tricky and was given incorrectly in The correct resu 
the total energy momentum density, derived first in j^], is, in the notation of Ji 



t, for 



P + A/2, S | I P - A/2, S } = AA bel P»P» 

-i^{[A bel + B bel ][P^™P + P»e 



u p/iafSi 



+ 



2A bel P»P" 
Pq + M 



.Opa/3 



J M' 



Pr 



(8) 



where A be i(A 2 ), B be i(A 2 ) are form factors, scalar functions of A 2 , and we have written 
A^i, B for A be i(0), B be i(0) respectively. Putting Eq. (jSJ) into Eq. (0), and taking P M = 
(Po, 0, 0, P), we obtain 
M<P,S| j£|P,S) 



(2tt) 3 5 3 (0) 



{(Aba + B be i)[2P 2 e jka0 - P(2P e jka3 + pV 0q3 - pi e k0a3 )} 



2AbelPPo rpk £ 0ja3 _ pj 0fcQ3l \ 

* M 



(9) 



P + M 

Since the separate quark and gluon momenta are not conserved, the analogue of Eq. ([8]) for 
them contains an additional term proportional to g^ v . Expressions only valid when such 
terms do not appear will be indicated by u (g liU = 0)"@. For the quark part of Tj^ one has 

( P + A/2, S | TZ ,\P- A/2, S ) = 4Am,, P"P" 



"4^ 



+ 



2A 



(<T = 0) 



Jdpafi j °Q p 



/9 



so that analogous to Eq. 
(2vr) 3 5 3 (0) 



P + M M 
the quark Belinfante angular momentum is given by 



(10) 



(A be i, q + B behq )[2P^ ka0 - P(2P e^ 3 + PV 0q3 - pJ e k0a3 )} 

(11) 



2^4&ei, qPPo j-pfc^0ja3 _ pj e 0ka3 



Po + M 



1^ 
J M' 



The connection with the notation in 



is : y4 = D/2,B = §-D/2 



2 Note that if light-cone Lorentz coordinates are permitted, this is not equivalent to "(fj, ^ v)" , since, for 
example g ++ = and = 1. 
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The operator form for the quark part Tj^ is, 



Let us define 



so that 



TZ, q = -jilvf D ^ + (/i i/)] (<T = 0) (12) 



If = JW^i (13) 



T g 9 = I [2^ + 37] or = o) (14) 



Because Tf is not symmetric under (/i <H- v) its matrix elements include an antisymmetric 
term, with an additional form factor 2t g = 2l g (A 2 = 0), and bearing in mind Eq. (I14j) . have 
the following structure: 



( P + A/2, <S | Sf | P - A/2, <S } = 4A 6e/ , P^P* 

+ B 6eU ] [P^ a/3 + P^e 



u p/ia/S] 



9 4 P 17 -\ *7 



9r ^ J ' P + M ^ J M' 

The operator form for the quark part of the Canonical Tj£L is 



can, q 

i - 
-i 



TZ hq =n4vfd v ik- (iT = 0) (16) 



and being non-symmetric under (/i -H- z/) its matrix elements have the form 
( P + A/2, S | | P - A/2, 5 ) = 4A cantq P"P v 



+ 2t can , [PV«* - PV^] + ri £ °^It7^ isT = 0) (17) 

rj) + 1V1 J IVl 



Now from Eqs. f ll2| 116]) . since (7 = 0, we have 



T£ q = d + M T^ q = -J n + d (is) 



Choosing the gauge A + = makes the two terms in Eq. (|18[) equal. Then, since is 
gauge invariant one has 

Tbel,q = ^6eZ,glgaugeA+=0 = ^am, q I gauge A+=0 (19) 

and thus, from Eqs. ( TTTJT) and ( |T71) . it is seen that 

Abel,q A-can, q | gaugeA+=0 and Bfyel^q P can, q\ gauge A+=0- (^0) 



Hence Eq. (fT5|) can be written 

( P + A/2, S | | P - A/2, S ) = 4A can , g | gaugeA+=0 P»P U 

■ f [ A I _i_ D | 1 rpA t , r P 1/Q /9 i p^PA 1 "^! 

4 |^ [-^mn, q \ gaugej4+=0 < J - > can, q |gauge^4+=0 J L ' J 

I Of Wl* cPva.fi — p^ c pua/3] I z '^can,g|gaugeA+=0 - 1 - 1 ^.Opaff I a p foi^ 

+ 4,^e ^ e j+ ^ + M e | M ^. ^i; 

Next, again choosing the gauge A + = and noting that is gauge invariant, we have 



2£ + = ^ + |gaugcA + =0 = ^lY D + ^| g augcA + =0 = Tg+ q \ smgeA+ =o (<? M+ = 0) (22) 
so that 2l g = ^can,q and thus the important consequence that when g^ u = 

( P + A/2, S | V? I P - A/2, S ) = ( P + A/2, 5 | T^, | P - A/2, 5 )| gaugeA+=0 . (23) 
Now, as shown by Shore and White [h| the following operator relation holds when g^ u = 
% X T bel q ~ % V T bel q = \ $-fiD v ^ - x v &fiD x i>] + - e» Xv(T ^ 7fT75 ij> + EoM + Divergence 



x A Xf - x u %f + - e^ va ^ 7ff 7 5 ^ + EoM + Divergence. (24) 



Ignoring the EoM terms which vanish via the equations of motion and the Divergence terms 
which vanish when taking matrix elements between localized states, and using Eq. (1231) 
yields the key relation 

( P + A/2, S | T b °l q - x k T°i l>g \ P - A/2, S ) 

= ( P + A/2, S | x>TZ, q - s*!^, | P - A/2, 5 )| gaugeA+=0 

-( P + A/2, 5 | ~ e ^ ^75 ^|P- A/2, 5 } (25) 

The LHS of Eq. (|25|) is just the matrix element of q (x) and the RHS can be recognized 
as the matrix element of 



M c l fc g (orbital;x) + M c l fc g (spin;x) = M^Jx) (gauged = 0) (26) 

so that 

M°£ q (x) = M^ fc g (x)| gaugeA+=0 . (27) 

Carrying out the Fourier transform and limiting process as in Eqs. E]) yields the remark- 
able result that in a spin 1/2 nucleon 

(P, S | JlX q \P,S)= (P, 5 I J^ q |P, 5 > | gauge A+=0- (28) 
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Now the total, quark plus gluon, Belinfante and Canonical angular momentum densities, as 
mentioned earlier, differ by a spatial divergence, so the total angular momenta 

pk _ pk pk , T jk . . T jk , T jk / 9q \ 

J bel ~ J bel,q ' J bel,G dJ - LU - J can ~ J can, q ' J can, G \ Z ' J ) 

may be taken as equal. It follows then from Eq. (128]) that also for gluons 

(P, S | j£ G \P,S) = (P, S I J^ G IP S ) | gaug M +=0 . (30) 

This quite unexpected result has interesting consequences for the several angular momentum 
sum rules in the literature 11], some based on the Belinfante, others on the Canonical 
versions of J, and this will be discussed in a paper in preparation. It also raises the intriguing 
question as to whether an analogous relation holds for a spin one hadron like the deuteron. 
It might be that the result only holds for the spin 1/2 case because of the limited structure 
of the matrix elements in e.g. Eq. (II ip . 
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